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ABSTRACT r—TA.  reactive-heat-pipe  concept  is  described  which  involves  the 
reaction  between  lithium  and  sulfur-hexafluoride.  For  this  concept,  the 
lithium  is  evaporated  from  the  inner  surface  of  the  wick  and  reacts  with 
gaseous  sulfur-hexafluoride;  the  energy  of  reaction  is  conducted  through 
the  wick  causing  lithium  to  evaporate  from  the  outer  surface  of  the  wick, 
the  vapor  then  transfers  the  energy  of  reaction  by  condensing  on  the  heat 
exchanger  surfaces  of  the  load.  A  small  combustor,  generating  7.5  kW, 
operating  at  condenser  temperatures  of  1200  K,  was  designed  and  tested  to 
demonstrate  the  operation  of  the  concept.  The  system  was  also  analyzed  in 
order  to  examine  the  design  characteristics  of  the  wick,  the  vapor  transport 
system  and  the  condensation  process  on  the  load.  The  results  indicate  an 
optimum  mesh  size  for  the  wick,  and  suggest  that  arterial  wicks  are 
desirable  for  some  operating  conditions.  The  greatest  thermal  resistance 
found  in  the  study  results  from  the  presence  of  noncondensibles  in  the  system; 
and  several  methods  of  control  are  suggested  to  minimize  this  difficulty. 
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If  hydrocarbon  fuels  are  used  to  fuel  the  power  source  of 
submersible  vehicles  for  deep-depth  operation,  there  are  problems 
encountered  due  to  volume  and  weight  requirements  and  the  inevitable 
complexity  of  design  and  construction  of  the  exhaust  system.  To 
avoid  these  difficulties  completely  closed  power  systems,  which  use 
thermal  energy  sources ,  have  been  developed  in  recent  years .  A 
thermal  energy  source  is  a  source  of  heat  at  an  elevated  temperature 
and  can  be  provided  by  exothermic  nuclear  or  chemical  reactions. 

The  energy  released  by  these  systems  can  be  used  to  supply  heat  for 
various  thermodynamic  power  cycles,  e.g.,  Brayton,  Stirling,  Rankine, 
etc. ,  which  ultimately  provide  mechanical  or  electrical  power. 

An  attractive  thermal  energy  source  for  submersibles  is  the 
liquid  metal  combustor.  The  heat  is  generated  by  the  chemical 
reaction  between  an  alkali  metal  and  a  gaseous  halogenated  oxidizer. 

In  the  normal  range  of  operating  temperatures,  the  reaction  products 
are  in  the  liquid  or  solid  state.  Since  the  volume  of  the  metal  salts 
produced  by  reaction  is  nearly  the  same  as  that  of  the  original  fuel, 
there  is  no  need  to  exhaust  the  products,  eliminating  the  problem  of 
product  disposal  while  underwater.  In  addition,  the  liquid  metal- 
halogen  gas  reactant  combinations  have  a  very  high  energy  density 
providing  a  compact  power  system. 


Pauluikonis  [1]  proposed  the  lithium-sulfur  hexafluoride 
reactant  combination  for  underwater  applications.  A  survey  of  many 
reactant  combinations  by  Van  der  Sluijs  [2]  also  leads  to  the 
conclusion  that  lithium-sulfur  hexafluoride  was  one  of  the  best  systems 
for  underwater  use.  These  reactants  have  extremely  high  reaction 
energy  output,  approximately  the  same  volume  of  fuel  and  products 
[3],  and  the  oxidant  is  nontoxic. 

The  reaction  is  completed  in  a  combustion  chamber  containing  the 
molten  fuel.  The  oxidizer  enters  the  chamber  as  a  gas  through  an 
injector.  The  stoichiometric  equation  for  the  lithium-sulfur 
hexafluoride  reaction  is 

8  Li (£)  +  SFfi(g)  +  6  LiF(Jl)  +  Li^U) 

The  products  of  the  reaction  are  liquids  for  the  usual  range  of 
operating  temperature  [3],  i.e.,  from  1100°  K  to  1200°  K.  The  average 
theoretical  energy  release  for  the  reaction  near  1200°  K  is  4.75 
kW-hr  per  kilogram  of  oxidizer  [3],  The  phase  equilibrium  between 
liquid  lithium  and  its  salts  yields  the  formation  of  an  immiscible 
liquid  region  and  the  molten  bath  splits  into  two  phases.  The 
product-rich  phase  settles  to  the  bottom  of  the  bath  due  to  its 
greater  density  and  can  be  accumulated  and  separated  easily. 

Thermal  energy  from  the  reaction  is  transferred  through  the  walls 
of  the  combustion  chamber  to  the  working  fluid  of  the  power  cycle, 
using  a  heat  exchanger  which  is  usually  integral  with  the  combustion 


chamber  walls. 


Uhlemann,  Spigt,  and  Hermans  [4],  Bierman  [5],  and  Faeth  and 
coworkers  [6-8]  have  investigated  the  metal  combustors  employing 
the  lithium-sulfur  hexafluoride  reactant  combination,  demonstrating 
many  of  the  operating  characteristics  of  these  systems.  The  overall 
objective  of  the  present  investigation,  is  to  study  a  novel  method  for 
conducting  this  reaction,  which  provides  an  improved  means  of 
transporting  the  thermal  energy  of  reaction  to  the  load. 

1,2  Previous  Related  Studies 

1.2.1  Submerged  Injector  Combustor.  One  method  of  conducting 
the  reaction  involves  the  use  of  a  submerged  injector  combustor  which 
is  illustrated  schematically  in  Figure  1.  The  exit  of  the  injector 
is  completely  submerged  in  the  liquid  fuel  bath.  The  oxidizer 
passes  into  the  combustion  chamber  through  this  injector.  The 
reaction  proceeds  in  a  turbulent  two-phase  jet  near  the  injector 
exit  [9].  The  gas  is  consumed  and  the  heavy  reaction  products 
settle  through  the  bath  and  are  collected  at  the  bottom  of  the 
chamber.  The  reaction  continues,  controlled  by  the  flow  rate  of 
oxidizer,  until  all  the  fuel  is  consumed.  At  this  point,  the  level 
of  the  product  layer  in  the  combustor  is  approximately  the  same  as 
the  initial  fuel  level  [3], 

There  are  several  disadvantages  introduced  by  direct  oxidizer 
injection  into  the  bath,  in  this  manner,  as  follows: 

i.  The  injector  is  subject  to  clogging  when  the  oxidizer 
flow  rate  is  varied  [7,  8], 
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Figure  1  Submerged  Injector  Combustor 


ii.  Due  to  the  proximity  of  the  molten  bath,  it  is  difficult  to 
cool  the  injector,  which  leads  to  corrosion  by  the  oxidizer, 
iii.  It  is  difficult  to  restart  the  combustor  since  the  injector 
clogs  shut  when  the  oxidizer  flow  is  terminated. 

Another  approach  which  can  resolve  some  of  the  problems  of  the 
submerged  injector  combustor  involves  the  coaxial  injector  combustor 
concept,  illustrated  schematically  in  Figure  2.  In  this  case,  the 
injector  has  two  passages,  one  for  a  core  flow  of  oxidizer  and  a 
second  outer  shroud  flow  of  inert  gas,  usually  argon.  The  inert  gas 
is  recirculated  through  a  vent  system.  The  outer  shroud  flow 
prevents  the  direct  contact  between  the  injector  body,  the  oxidizer, 
and  the  fuel  that  occurs  in  single  passage  injectors.  This 
eliminates  the  possibility  of  unstable  slugging  flow  which  occurs 
when  a  reacting  gas  is  injected  into  a  liquid  at  a  low  flow  rate. 

This  type  of  flow  instability  provides  the  mechanism  for  plugging 
the  flow  passage  with  product  material  at  low  flow  rates  for  a 
single  passage  injector  [8].  However,  problems  with  corrosion  of 
the  injector  still  exist.  In  addition,  the  coaxial  injector  system 
is  more  complicated  than  the  single  passage  injector,  tending  to 
reduce  reliability. 

Another  problem  encountered  with  submerged  injector  systems 
involves  the  variation  in  bath  heat  transfer  characteristics  as  the 
product  accumulates  in  the  combustor.  Unless  stirred,  the  product 
acts  as  a  stable  stratified  layer,  and  the  poor  heat  transfer 
characteristics  of  this  region  limits  capabilities  for  transferring 


heat  to  the  load. 
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1.2.2  Wick-Type  Combustor.  The  coaxial  injector  meets  many 
objectives  for  long-term  variable  load  operation;  however,  the  need 
to  vent,  filter,  and  recirculate  the  argon  flow  to  the  injector 
increases  system  complexity.  In  order  to  achieve  the  full 
capabilities  of  the  lithium/sulf ur  hexafluoride  reaction  concept, 
a  system  which  requires  no  argon  flow  must  be  considered.  This  can 
be  done  with  the  wick-type  combustor,  which  is  illustrated 
schematically  in  Figure  3  [10].  The  arrangement  is  such  that  the 
lithium  containing  wick  is  mounted  directly  on  the  inside  surface  of 
the  combustor,  and  extends  below  the  level  of  the  bath.  The  wick  can 
control  liquid  motion  near  the  injector,  in  order  to  circumvent 
liquid  slugging  in  the  injector  passage  which  led  to  clogging  of 
single-passage  injectors  at  low  flow  rates  [6],  Capillary  forces 
draw  lithium  through  the  wick  from  the  bath.  The  oxidizer  passes 
through  an  injector  above  the  liquid  surface  into  the  ullage  region. 
The  reaction  proceeds  as  the  lithium  evaporates  from  the  wick  and 
the  energy  of  reaction  is  transferred  to  the  load  through  the  wicking. 
The  products  condense  and  fall  into  the  liquid  bath  at  the  bottom  of 
the  combustor. 

Blakeslee  [10]  demonstrated  that  the  combustion  process  occurs 
in  the  gas  phase  near  the  surface  of  the  wick,  as  pictured  in  Figure 
4.  Fuel  evaporates  from  the  wick  and  diffuses  toward  the  reaction 
zone,  while  the  oxidizer  diffuses  from  the  ullage  space  to  the 
reaction  zone.  The  reaction  proceeds  as  a  diffusion  flame,  however, 
the  products  condense  as  they  diffuse  to  cooler  regions  outside  the 
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Figure  3  Wick-Type  Combustor 
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reaction  zone.  Some  products  condense  on  the  surface  of  the  wick,  as 
well,  in  a  dropwise  manner  [10]. 

For  combustion  rates  observed  during  wick-type  burning,  gas 
velocities  normal  to  the  surface  of  the  wick  are  small  and  the 
pressure  of  the  process  is  essentially  constant.  This  implies  that 
the  pressure  within  the  ullage  region  is  fixed  by  the  vapor  pressure 
of  the  fuel  at  the  surface  of  the  wick.  Figure  5,  taken  from 
Reference  [10],  is  an  illustration  of  the  comparison  between 
measured  pressures  within  the  combustor  and  predictions  based  on  the 
vapor  pressure  of  the  fuel  at  the  surface  temperature  of  the  wick. 
Except  at  the  lowest  pressures  where  effects  of  noncondensibles 
and  experimental  errors  are  greatest,  the  agreement  between  the  vapor 
pressure  prediction  and  the  measured  pressure  is  reasonably  good. 

Heat  transfer  from  the  reaction  zone  to  the  load  occurs  through 
the  wick,  and  is  not  influenced  by  the  amount  of  product  within  the 
combustor.  Injector  reliability  is  also  improved  over  submerged 
injector  designs  [4],  and  the  complexity  of  a  recirculating  argon 
flow  is  removed.  While  these  features  are  improvements  over  the 
submerged  injector,  several  problems  remain,  as  follows: 

i.  Since  the  wick  is  fixed  to  the  inside  surface  of  the 
combustion  chamber,  the  wick  is  difficult  to  replace, 

ii.  The  heat  transfer  characteristics  of  the  load  strongly 

influence  wick  temperatures,  excessive  local  cooling  can 
cause  the  wick  to  be  progressivley  extinguished, 

iii.  The  configuration  is  relatively  inflexible  since  the 

wick  must  cover  the  interior  surfaces  ir.  ^de  to  prevent 
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corrosion.  In  order  to  avoid  wick  assembly  problems 


practical  designs  tends  to  use  excessivle  ullage  volumes, 
reducing  system  compactness. 

1.3  Reactive-Heat-Pipe  Concept 

In  order  to  overcome  the  problems  in  the  wall-mounted-wick - 
combustor,  a  reactive-heat-pipe-combustor  is  proposed.  This  system 
combines  the  features  of  the  wick-type-combustor;  and  a  conventional 
heat  pipe,  for  heat  transport  capabilities.  The  reactive-heat-pipe- 
combustor  is  illustrated  in  Figures  6  and  7.  The  heat  exchanger  is 
located  outside  the  combustor  in  Figure  6  and  inside  the  combustor 
in  Figure  7. 

The  wick  is  free-standing  in  the  reaction  vessel,  attached  at 
its  upper  end  at  the  top  of  the  combustion  chamber,  and  opening  into 
the  bath  at  the  bottom.  The  oxidizer  is  injected  into  the  space 
enclosed  by  the  wick  and  bath.  Reaction  is  confined  within  this 
volume,  and  proceeds  in  the  same  manner  as  the  wall-mounted  wick 
arrangement. 

The  energy  of  reaction  is  transported  to  the  inside  surface  of 
the  wick  by  convection,  radiation,  and  condensation  of  the  product. 
Then  the  heat  is  conducted  through  the  wick  and  causes  lithium  to 
evaporate  from  the  outside  surface  of  the  wick.  The  vapor  leaving 
the  wick  flows  to  the  heat  exchanger  through  the  gas  space  and  then 
condenses  to  complete  the  transfer  of  reaction  energy.  The  condensate 
is  returned  to  the  molten  storage  bath  in  the  combustion  chamber. 

Since  the  temperature  of  the  inside  surface  of  the  wick  is  higher  than 
the  outside  surface,  there  is  a  pressure  difference  across  the  wick. 


Lth  E: 
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This  pressure  difference  is  accommodated  by  a  slight  displacement  of 
the  liquid  level  under  the  reaction  zone. 


There  are  several  advantages  for  this  reactive-heat-pipe 
configuration  as  follows: 

i.  The  heat  pipe  for  energy  transfer  to  the  load  helps  to 

maintain  relatively  uniform  temperature  distribution  on  both 
the  wick  and  the  load  heat  exchanger.  This  simplifies 
the  design  of  the  heat  exchanger  by  minimizing  thermal 
stresses . 

ii.  Since  the  wick  is  free-standing  in  the  combustion  chamber, 
its  replacement  is  convenient. 

iii.  The  wick  is  independent  of  the  load,  therefore,  danger  of 
local  subcooling  and  extinguishment  is  eliminated, 

iv.  The  wick  confines  the  reaction  zone  to  a  small  volume 

within  the  combustion  chamber,  allowing  the  heat  exchanger 
to  be  located  in  the  ullage  space.  This  conserves  system 
volume  and  improves  compactness, 

v.  The  heat  pipe  action  allows  the  energy  of  reaction  to  be 

transported  some  distance  from  the  reaction  area  improving 
the  flexibility  of  the  system  layout, 

vi.  The  wick  and  the  load  heat  exchanger  have  relatively  small 
heat  capacities  providing  rapid  adjustment  to  load  changes 
and  a  fast  start.  The  ignitor  system  can  also  be  located 
near  the  wick  so  that  just  a  portion  of  the  wick  can  be 


ignited. 


1.4  Specific  Objectives  of  the  Stud 


As  discussed  above,  the  reactive-heat-pipe-combustor  has  several 
advantages  over  earlier  liquid  metal  combustors.  It  is  capable  of 
fast  reaction  to  load  changes;  relatively  fast  start;  flexible  heat 
exchanger  configuration;  independent  operation  of  the  wick  and  the 
load  heat  exchanger;  capability  for  complete  shutdown  and  restart  any 
number  of  times;  and  uniform  temperature  distributions  on  both  the 
wick  and  the  load,  which  simplifies  the  design  of  both  these 
components . 

While  earlier  work  [10]  has  suggested  the  feasibility  of  the 
reactive-heat-pipe  concept,  this  configuration  has  not  yet  been 
operated  and  its  design  characteristics  have  not  been  studied. 
Therefore,  the  specific  objectives  of  the  present  study,  are  as 
follows: 

i.  Design  and  test  a  small  scale  version  of  the  reactive-heat- 
pipe-combustor  in  order  to  observe  its  practical  operation, 
ii.  Complete  a  theoretical  analysis  of  this  type  of  combustion 
system.  The  purpose  of  the  analysis  will  be  to  obtain 
design  characteristics  of  the  wick  and  the  heat  transport 
system,  and  to  investigate  their  performance  characteristics 

The  study  is  limited  to  the  lithium-sulfur  hexafluoride 


reactant  combination.  The  principles  can  be  applied  to  other 
reacting  combinations  in  this  family,  however,  these  reactants 
appear  to  have  the  greatest  advantage  for  practical  applications. 


CHAPTER  II 


EXPERIMENTAL  APPARATUS 

2.1  General  Characteristics 

A  combustor  was  designed,  fabricated,  and  tested  in  order  to 
determine  the  operating  characteristics  of  the  reactive-heat-pipe 
configuration.  In  order  to  conserve  costs  the  combustor  was 
relatively  small,  containing  about  .9  kg  of  lithium  and  generating 
a  maximum  thermal  power  of  7.5  kW.  The  heat  exchanger  for  the 
combustor  was  also  an  inexpensive  radiative-type  arrangement.  The 
load  provided  by  the  heat  exchanger  was  controlled  by  moving  a 
radiation  shield  over  the  surface  of  the  heat  exchanger. 

2.2  Combustor 

A  sketch  of  the  reactive-heat-pipe-combustor  is  shown 
schematically  in  Figure  8,  the  major  features  of  the  design  are 
summarized  in  Table  1.  The  combustor  body  is  cylindrical  in  shape, 
constructed  of  AISI  Type  316  stainless  steel.  The  cylinder  is 
mounted  in  an  upright  position  and  supported  from  its  top  flange. 

The  combustor  is  filled  about  half-way  with  fuel,  the  lower  portion 
acting  as  a  fuel/product  storage  reservoir,  the  upper  portion  acting 
as  the  combustion  area  and  heat  exchanger. 

The  oxidizer  enters  the  combustor  through  a  nickel  injector 
located  at  the  center  of  the  top  flange  of  the  combustor.  The  flow 
passage  of  the  injector  had  a  diameter  of  10  mm.  The  flange  area  near 
the  exit  of  the  injector  is  reduced  in  thickness  in  order  to  minimize 


Figure  8  Reactive-Heat-Pipe  Combustor 
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the  heat  flow  into  this  region  and  help  to  keep  the  exit  area  cool 
in  order  to  avoid  corrosion  of  the  metal  by  SF^  at  high  temperatures . 
The  injector  also  incorporated  a  widened  cylindrical  portion  which 
provided  natural  convection  cooling.  Nickel  is  used  in  this 
component  in  order  to  reduce  the  corrosive  action  of  the  oxidizer 
at  high  temperatures . 

The  combustor  was  sealed  by  welding  the  outer  edge  of  the  top 
flange  assembly.  This  edge  was  extended  so  that  it  could  be  ground 
away,  allowing  the  combustor  to  be  opened  for  inspection  and 
refueling  several  times  before  replacement  was  necessary. 

The  wick  is  spot-welded  coaxially  to  the  top  flange  of  the 
combustor  and  hangs  down  into  the  bath.  The  wick  is  constructed  of 
two  layers  of  100-mesh  type  316  stainless  steel  screen.  The  lower 
end  of  the  wick  extends  to  a  position  about  50  mm  from  the  bottom  of 
the  combustor.  The  active  area  of  the  wick,  above  the  liquid 
surface,  has  a  length  of  approximately  125  mm. 

The  cylindrical  combustor  wall  above  the  liquid  level  acts  as 
the  load  heat  exchanger.  Heat  is  transferred  from  this  surface 
primarily  by  radiation,  with  a  small  additional  contribution  from 
natural  convection  at  typical  operating  temperatures.  The  rate  of 
heat  transfer  from  the  heat  exchanger  is  controlled  by  varying  the 
vertical  position  of  an  adjustable  radiation  shield.  The  completely 
closed  position  of  the  shield  corresponds  to  the  no-load,  idling 
condition. 

During  the  course  of  testing  it  became  evident  that  a  baffle  over 
the  inside  condensing  surface  of  the  heat  exchanger  would  be 


desirable  in  order  to  control  the  heat  transfer  resistance  due  to 
noncondensible  gases  in  the  system.  The  baffle  is  attached  to  the 
top  of  the  combustor  and  passes  around  the  wick  in  a  coaxial  manner. 
The  flow  of  vapor  from  the  wick  to  the  condensing  surface  of  the 
heat  exchanger  is  illustrated  in  Figure  9  for  both  the  unbaffled 
and  baffled  arrangements.  The  action  of  the  baffle  is  to  force  the 
vapor  flow  into  a  tangential  path  along  the  surface  of  the  heat 
exchanger  sweeping  the  noncondensible  gas  toward  a  vertical  collection 
tube  located  in  the  top  flange  of  the  combustor.  The  noncondensible 
gas  can  be  vented  from  this  gas  collector  in  a  manner  that  will  be 
described  later. 

The  combustor  could  be  heated  electrically  in  order  to  melt  the 
fuel  and  start  the  combustion  process.  The  electrical  start  heaters 
are  located  in  the  movable  radiation  shield.  They  consist  of  two 
bead  heaters  manufactured  by  the  Cole-Parmer  Instrument  Company 
(Model  316-40) ,  providing  an  output  of  150  watts  per  meter  of 
length.  The  total  maximum  power  of  the  two  heaters  was  2.7  kW.  The 
heating  power  was  regulated  by  two  VARIC  auto  transformers.  The 
electrical  heaters  are  insulated  externally  by  the  insulation  in  the 
radiation  shield. 

2.3  Gas  Supply  System 

A  sketch  of  the  gas  supply  system  is  shown  in  Figure  10.  The 
entire  system  is  welded  or  brazed  wherever  possible  to  reduce  the 
possibility  of  leakage.  Those  connections  which  could  not  be  welded 
or  braised  were  sealed  with  a  silicon-rubber  sealant.  A  VEECO  MS-90AB 


CONDENSER  WALL 


Figure  9  Top  View  of  Condenser  Baffle  Arrangement 
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Figure  10  Schematic  of  Gas  Supply  System  for  Reactive-Heat-Pipe  Combustor  Tests 
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leak  detector  was  used  to  helium-leak  check  the  entire  system  before 
each  test. 

The  oxidizer  is  stored  as  a  liquid  in  a  cylinder  and  drawn  off 
from  the  gas  phase  through  a  pressure  regulator.  The  oxidizer 
passes  through  a  critical  flow  orifice  which  is  used  to  meter  the 
flow  and  then  through  the  injector  into  the  reaction  region.  A  line 
near  the  inlet  of  the  injector  goes  back  to  a  mercury  manometer  and 
a  thermocouple  gage  in  order  to  record  the  gas  pressure.  The  line 
to  the  pressure  gages  contains  a  trap  which  is  used  to  protect  the 
gages  from  contact  with  lithium  vapors. 

The  combustor  is  connected  to  a  vacuum  pump  through  valve  4  on 
the  condenser  side  of  the  wick  and  valve  2  on  the  reaction  side  of 
the  wick.  These  connections  allow  noncondensibles  to  be  removed  from 
the  system.  The  helium-leak  detector  is  also  connected  to  the  vacuum 
side  of  the  system  through  valve  3  and  the  load  heat  exchanger. 

2.4  Instrumentation 

The  oxidizer  flow  to  the  combustor  was  controlled  by  a  pressure 
regulator  and  measured  with  a  calibrated  critical-flow  orifice.  The 
pressure  regulator  was  a  Matheson  Model  9  two-stage  regulator  with 
10-250  psig  capability.  The  upstream  pressure  of  the  critical-flow 
orifice  was  measured  with  a  Heise  absolute  pressure  gauge  (0-300  psia 
range  with  0.5  psi  divisions).  The  orifice  was  calibrated  with  a 
Precision  Scientific  Wet-Test  Meter  (0.1  cubic  feet/revolution). 

Combustor  pressures  were  measured  with  a  thermocouple  vacuum 
gauge,  Model  ML-25,  made  by  Mill  Lane,  Engineering  Co.,  Inc.,  and  an 
absolute  mercury  manometer  with  1  mm  subdivisions,  made  by  Meriam 


Instrumentation  Co.  The  pressures  were  measured  at  the  inlet  of  the 
injector.  Since  the  oxidizer  flows  used  during  tests  were  relatively 
low  (less  than  1.5  kg/hr)  the  pressure  drop  through  the  injector  was 
negligible  [10].  Therefore,  within  the  accuracy  of  the 
instrumentation,  the  measured  pressure  was  equal  to  the  combustor 
pressure . 

The  surface  temperature  of  various  components  in  the  apparatus 
were  measured  with  thermocouples .  The  thermocouples  were 
constructed  with  chromel-alumel  wire  manufactured  by  the  Thermo- 
Electric  Company,  type  G/G-26-K,  within  American  National  Standards 
Institute  limits  of  error  of  +  .75  percent  in  the  temperature  range 
of  use.  The  thermocouples  were  spot  welded  directly  to  the 
combustor  surfaces.  There  were  twenty  thermocouples  welded  to  the 
combustor  wall,  6  to  the  injector,  6  to  the  surface  of  the  upper 
cover,  and  4  to  the  noncondensibles  trap.  The  thermocouples  on  the 
combustor  wall  were  located  in  three  lines,  equally  spaced  around  the 
periphery  of  the  combustor , with  one  line  directly  below  the  inlet 
of  the  noncondensible  gas  trap.  The  vertical  spacing  of  the 
thermocouples  was  26  mm  in  the  vapor  region,  and  78  mm  in  the  liquid 
region. 

Temperatures  were  recorded  every  eighteen  seconds  on  three 
Leeds  and  Northrup  Speedomax-H  twelve-point  recorders.  The  extension 
leads  and  recorders  were  calibrated  with  a  Leeds  and  Northrup  model 
8686  millivolt  potentiometer.  The  temperature-millivolt  conversion 
values  were  taken  from  National  Bureau  of  Standards  tables. 


2.5  Reactants 


The  lithium  used  for  the  tests  was  obtained  from  Lithium 
Corporation  of  America.  The  fuel  was  packed  in  cans  and  sealed  under 
argon.  The  lithium  is  guaranteed  as  99.9  percent  pure  with  the  major 
impurity  being  sodium. 

The  sulfur  hexafluoride  was  obtained  from  Matheson  Gas  Products. 
The  gas  is  contained  in  a  cylinder  as  a  liquid  under  its  own  vapor 
pressure  (2200  KPa  at  294  K) .  The  gas  was  guaranteed  to  have 
minimum  purity  of  99.8  percent. 

2.6  Experimental  Procedure 

Prior  to  loading  the  fuel,  the  interior  of  the  combustion 
chamber  and  the  upper  cover,  which  is  integral  with  injector  and 
the  noncondensibles  collector,  were  rinsed  thoroughly  with  methanol, 
then  dried  in  a  vacuum/glove  box.  The  fuel  loading  steps  are  as 
follows:  (1)  the  vacuum  box  was  filled  with  argon  until  it  reached 

atmospheric  pressure;  (2)  the  lithium  was  cleaned  and  weighed  in 
the  box  to  yield  a  fuel  loading  of  around  0.9  kg;  (3)  the  combustor 
was  heated  to  525  K  in  the  vacuum  box,  melting  the  lithium  in  the 
combustion  chamber;  (4)  after  melting  the  lithium,  the  upper  cover 
was  aligned  evenly  with  the  combustor  body  and  the  flanges  of  the 
cover  and  the  body  were  clamped  together;  (5)  the  loaded  combustor 
was  then  cooled  to  room  temperature  and  the  venting  holes,  and  the 
seam  between  the  flanges  of  the  cover  and  the  body  were  sealed  with 
temporary  seal;  (6)  the  combustor  was  then  removed  from  the  vacuum 
box  and  the  top  cover  was  seal-welded  to  the  combustor  body.  During 


welding  a  continuous  argon  purge  was  maintained  into  the  combustor 
to  avoid  contaminating  the  lithium  with  air. 

After  welding,  the  combustor  was  connected  to  the  oxidizer  supply 
line  and  the  vent  line.  Then  the  system  was  checked  for  leaks. 

A  test  was  begun  by  heating  the  combustor  under  continuous  vacuum 
to  about  850°  K  to  degas  the  fuel.  At  this  point,  the  vacuum  pump 
was  shut  off,  but  heating  continued.  When  the  temperature  of  the 
bath  reached  1100°  K,  a  small  amount  of  oxidizer  was  bled  in  the 
combustor.  When  chemical  reaction  began,  the  pressure  of  the 
combustor  would  suddenly  drop  and  a  sharp  temperature  rise  of  the 
combustor  wall  occurred.  After  ignition  the  oxidizer  was  supplied 
steadily  and  the  heater  was  shut  off. 

Oxidizer  flow,  combustor  pressure,  and  system  temperature  were 
recorded  from  the  ignition  point  to  system  shut  down,  at  intervals 
of  1,  2,  or  5  minutes.  The  combustor  was  vented  to  vacuum  when 
excessive  noncondensibles  accumulated  in  the  system. 

Operation  was  ended  by  gradually  terminating  the  oxidizer  flow 
and  allowing  the  combustor  to  cool  to  room  temperature.  In  order  to 
prevent  reaction  with  air  in  the  event  of  possible  leakage,  the 
combustor  ullage  space  was  filled  with  argon  when  the  combustor  was 


stored. 


CHAPTER  III 


EXPERIMENTAL  RESULTS  AND  DISCUSSION 

Six  reactive-heat-pipe  combustor  tests  were  conducted,  Table  2 
summarizes  the  overall  characteristics  of  these  tests.  The  main 
objectives  of  the  tests  were  to  demonstrate  the  feasibility  of  the 
concept,  to  check  the  operation  of  the  combustor,  and  to  improve  the 
design  on  the  basis  of  observed  operating  characteristics.  In  the 
following  the  results  for  each  test  are  discussed  in  detail. 

3.1  Test  1 

This  test  was  conducted  to  gain  familiarity  with  the  test 
apparatus  and  the  combustor.  No  attempt  was  made  to  achieve  good 
fuel  utilization,  therefore,  the  fuel  utilization  in  this  test  is 
rather  low.  The  combustor  was  operated  for  a  period  of  4.3  hours, 
and  then  shut  down  as  planned.  No  difficulties  were  encountered  with 
respect  to  injector  operation  or  clogging.  The  major  problem  during 
the  test  was  insufficient  condenser  action  which  prevented  operation 
at  high  heat  flux.  Therefore,  the  power  density  on  the  wick  for  this 
test  is  rather  low.  The  poor  condenser  action  was  attributed  to  the 
presence  of  small  quantities  of  noncondensible  gas  between  the  wick 
and  the  condenser  surface.  The  quantities  were  small,  since  combustor 
pressure  levels  were  normal  for  the  temperature  levels  of  the 
combustor.  Even  small  quantities,  however,  can  inhibit  condenser 
action,  since  the  vapor  sweeps  the  noncondensible  gas  against  the 


Run  Wick  Power  Wall  Fuel 

Test  Length  Density  Temperature  Utilization  Type 

(hours)  (W/cm2)  (°  K)  (%) 


condenser  surface  providing  a  diffusional  barrier  to  the  condensation 
process . 

The  combustor  was  disassembled  after  the  test  in  order  to 
examine  the  interior  components.  No  corrosion  was  evident  in  the 
injector  area.  The  wick  was  wetted  (to  the  top)  with  lithium 
verifying  the  design  of  the  wick  assembly.  The  wick  was  also  in 
good  condition,  with  no  evidence  of  corrosion. 

These  results  verify  the  general  operation  of  the  reactive  heat 
pipe  combustor,  and  the  design  procedures  used  for  wick  fabrication. 
The  results  also  indicated,  however,  that  corrective  action  was 
required  in  order  to  reduce  the  heat  transfer  resistance  of  the 
condenser  due  to  the  presence  of  noncondensibles. 

3.2  Test  2 

In  order  to  control  the  quantities  of  noncondensible  gas  within 
the  condenser  area,  a  vertical  condenser  tube  was  installed  for  this 
test,  as  shown  in  Figure  8.  A  vent  line  at  the  upper  end  of  the 
condenser  tube  allowed  noncondensible  gases  to  be  removed  from  the 
system. 

When  noncondensible  gas  was  present  in  the  condenser  region,  it 
impeded  heat  transfer  and  cool  dark  spots  appeared  on  the  combustor 
wall,  particularly  near  the  noncondensible  gas  trap.  After  the  gas 
was  vented  from  the  trap,  the  combustor  wall  became  uniformly 
bright  red.  It  was  only  necessary  to  remove  gas  at  the  start  of  the 
test;  no  additional  accumulation  was  observed  as  the  test  proceeded. 
The  amount  of  noncondensible  gas  in  the  trap  could  be  determined  by 


measuring  the  length  of  the  cool  dark  area  at  the  top  of  the 
condenser  tube,  where  condensing  action  was  inhibited. 

Figure  11  is  a  plot  of  the  oxidizer  flow  rate,  the  combustor 
pressure,  and  the  average  temperature  of  the  condenser  wall  during 
the  test.  The  oxidizer  flow  rate  was  adjusted  to  maintain  the 
condenser  wall  temperature  near  1150  K,  and  the  load  was  varied  by 
shifting  the  position  of  the  heat  shield.  Since  the  combustor 
temperature  was  nearly  constant,  the  combustor  pressure  did  not  vary 
greatly  as  the  load  was  changed. 

During  the  high  load  period,  beginning  at  180  minutes,  the 
condenser  wall  temperature  dropped  somewhat,  and  the  combustor 
pressure  increased,  indicating  an  increase  in  the  temperature  of  the 
wick.  Due  to  the  greater  heat  flux  through  the  wick  at  this 
condition,  some  increase  in  combustor  pressure  is  expected  at  high 
loads.  The  pressure  change  observed  in  Figure  11,  however,  is 
greater  than  is  predicted  from  Figure  5  for  these  operating 
conditions.  This  observation  suggests  that  small  quantities  of 
noncondensibles  were  still  present  on  the  main  condensing  surface  of 
the  combustor,  even  though  the  vertical  condenser  tube  removed  most 
of  the  gas. 

Figure  12  illustrates  the  temperature  distribution  along  the 
combustor  wall  at  various  load  conditions  (note  that  distance  is 
plotted  from  the  top  of  the  combustor,  so  that  the  condenser  area  is 
in  the  lower  portion  of  each  figure  but  above  the  liquid  level.  At 
full  load  conditions,  the  condenser  wall  temperature  is  nearly 
constant,  which  illustrates  the  desirable  characteristics  of  the 
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heat-pipe  principle  utilized  in  this  combustor.  Part  load  is 
achieved  by  removing  the  radiation  shield.  Under  these  conditions 
the  upper  portion  of  the  condenser  is  operating  at  full  load,  while 
the  lower  portion  is  essentially  at  an  idle  condition.  This  leads  to 
some  loss  of  temperature  uniformity,  although  the  variation  does  not 
exceed  20  K.  Even  this  variation  is  somewhat  larger  than  would  be 
predicted,  due  to  the  presence  of  small  amounts  of  noncondensible 
gas,  as  discussed  earlier. 

The  combustor  was  operated  until  pressures  began  to  rise, 
indicating  that  the  fuel  was  nearly  consumed.  At  this  point  the 
operation  of  the  combustor  was  terminated.  The  overall  fuel 
utilization  of  this  test  was  96%,  which  is  quite  satisfactory. 

After  disassembling  the  combustor,  inspection  of  the  interior 
components  indicated  that  the  injector  and  the  wick  were  in 
excellent  condition. 

3.3  Test  3 

The  combustor  arrangement  for  this  test  was  the  same  as  Test  2. 
This  test  was  somewhat  shorter  than  the  previous  test,  since  there 
were  two  one-half  hour  periods,  at  full  load,  as  opposed  to  one 
high  power  period  in  Test  2.  Operation  was  terminated  when  fuel 
was  consumed,  yielding  an  overall  fuel  utilization  of  97%. 

The  test  was  relatively  free  of  problems,  no  corrosion  or 
clogging  was  observed,  and  wall  temperatures  at  full  load  conditions 
were  uniform — similar  to  the  results  of  Test  2.  Actual  condenser 
performance,  however,  was  not  as  good  as  predicted,  due  to  the 
continued  presence  of  small  amounts  of  noncondensible  gas  in  the 
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system.  The  heat  transfer  coefficient  of  the  condenser  was  calculated 
as  described  in  Appendix  A.  It  was  found  that  the  heat  transfer 
coefficient  on  the  condenser  was  lower  than  the  expected  value  by  a 
factor  of  24.  This  strongly  suggested  that  the  design  of  the 
condenser  area  required  improvement  in  order  to  achieve  the  full 
potential  of  the  system. 

3.4  Test  4 

In  order  to  further  control  noncondensibles,  a  baffle  was 
installed  as  shown  in  Figure  9.  The  baffle  was  designed  to  prevent 
an  excessive  pressure  drop  in  the  condenser. 

When  a  baffle  was  not  present,  the  vapor  flow  sweeps 
noncondensibles  to  the  trapped  area  at  light  loads .  For  full  load 
conditions,  however,  the  primary  vapor  flow  directly  to  the  vapor 
surface  competes  with  the  flow  to  the  trap  and  noncondensibles  can 
collect  along  the  wall.  The  baffle  introduces  a  tangential  vapor 
flow  which  now  tends  to  sweep  the  wall  clear  of  noncondensibles  at 
full  load,  forcing  the  gases  toward  the  trap  where  they  do  not  impede 
condenser  action. 

The  baffle  worked  as  planned  during  the  early  portions  of  the 
test,  pushing  the  noncondensibles  to  a  narrow  region  near  the  trap. 
However,  due  to  an  excessively  light  construction,  the  baffle  broke 
due  to  thermal  cycling  and  failed  during  the  latter  portions  of  the 
test  (this  was  observed  as  spread  of  the  dark  region  uniformly 
around  the  combustor  and  from  the  condition  of  the  baffle  following 


the  test). 


In  most  other  aspects,  this  test  was  similar  to  Test  3. 

Reasonably  good  fuel  utilization  was  obtained,  and  the  power  density 
of  the  wick  could  be  varied  up  and  down  over  the  range  idle  -  26.0 
W/sq-cm.  The  operation  of  the  heat  pipe  was  not  affected  by  fuel 
consumption  until  very  near  the  end  of  the  test. 

While  the  vent  on  the  noncondensible  trap  could  be  used 
satisfactorily  in  the  early  portion  of  the  test,  the  line  tended  to 
become  clogged  with  lithium  as  the  test  proceeded.  This  behavior 
indicated  the  need  for  some  modification  of  the  trap,  if  continued 
removal  of  noncondensible  gases  was  desired  throughout  combustor 
operation. 

3.5  Tests  5  and  6 

These  tests  were  run  to  examine  capabilities  for  restarting  the 
combustor.  The  trap  was  also  modified  by  placing  it  in  a  horizontal 
position,  in  an  attempt  to  reduce  blockage  of  the  vent  line  by 
condensed  lithium,  which  was  pumped  to  the  top  of  the  tube  by 
natural  convection  instabilities  when  the  vertical  arrangement  was 
used. 

The  first  test.  Test  5,  was  run  for  a  period  of  three  hours, 
over  a  variable  load  cycle,  with  two  high-power  periods,  each  one-half 
hour  long,  and  then  shut  down  as  planned.  The  second  test.  Test  6, 
was  a  restart  of  the  first  test,  after  a  storage  period  of  three  days. 
The  combustor  was  operated  for  1.3  hours  during  the  second  test  and 
shut  down  upon  indication  of  consumption  of  fuel  due  to  increasing 
combustor  pressure. 


The  overall  fuel  utilization  for  the  two  tests  was  in  excess  of 
95%.  The  test  series  demonstrated  the  feasibility  of  restarting 
the  reactive-heat-pipe-combustor . 

Injector  operation  was  satisfactory  and  no  thermal  cycling 
problems  were  encountered  with  the  combustor  materials. 

Noncondensible  gas  could  be  removed  early  in  the  period  of  operation 
from  the  vent  in  the  gas  collector.  However,  the  small  vent  line 
eventually  clogged  with  condensed  lithium.  Fortunately  the  remaining 
levels  of  noncondensibles  were  low,  so  that  system  operation  could  be 
undertaken  at  full  power  levels  with  the  low  levels  of 
noncondensible  gas  remaining  in  the  system  stored  in  the  volume  of 
the  trap.  The  flow  of  oxidizer  did  not  introduce  additional 
quantities  of  noncondensibles,  in  detectable  amounts;  which  was 
true  for  all  the  tests.  The  interior  components,  wick,  baffle, 
injector,  etc.,  were  in  excellent  condition  following  the  test.  No 
difficulty  was  encountered  with  failure  of  the  baffle,  due  to  its 
heavier  construction  for  this  test. 

3.6  Discussion 

The  test  series  undertaken  during  this  investigation 
demonstrated  that  the  reactive-heat-pipe-combustor  operates  as 
anticipated.  The  major  problem  encountered  during  the  tests, 
involved  the  presence  of  noncondensibles  within  the  region  of  the 
condenser,  which  increased  the  heat  transfer  resistance.  From  the 
results  of  Test  5  and  Test  6,  the  amount  of  noncondensibles  can  be 
controlled  by  properly  degassing  the  fuel  and  the  combustor  walls; 
the  remaining  small  amounts  of  noncondensibles  can  be  controlled 


CHAPTER  IV 


THEORETICAL  MODEL 

Experimentation  with  a  combustor  having  a  specific 
configuration  is  of  value  for  demonstrating  the  operation  of  the 
concept.  However,  this  approach  is  limited  in  that  it  cannot 
consider  all  aspects  of  system  design  without  excessive  testing. 
Therefore,  in  the  following,  an  analysis  of  the  reactive-heat-pipe 
system  is  discussed,  in  order  to  study  the  general  operating 
limitations  of  the  concept  over  a  wider  range  of  variables  than  was 
accessible  experimentally.  In  many  ways,  the  system  is  similar  to 
a  heat  pipe,  and  the  present  analysis  has  made  liberal  use  of 
theoretical  findings  for  heat  pipes.  New  phenomena  that  must  be 
considered  involve  reaction  from  one  surface  of  the  wick  and  the 
effect  of  the  pressure  difference  across  the  wick. 

The  analysis  addresses  three  potential  failure  modes  as  follows: 

1.  The  wick  can  fail  if  it  is  unable  to  transfer  sufficient 

lithium,  to  the  top  of  the  wick,  for  both  reaction  and  heat 
transport  requirements.  Temperatures  within  the  wick  are 
maintained  at  acceptable  levels,  in  the  presence  of 
reaction,  by  lithium  evaporation  from  the  surfaces  of  the 
wick.  If  fuel  does  not  reach  the  upper  portions  of  the  wick, 
this  region  can  reach  high  temperature  levels  where  the 
material  is  directly  attacked  by  the  oxidizer,  resulting 
in  burn-out.  This  limit  is  referred  to  as  the  wick 


pumping  limit. 


2. 


The  system  will  fail  to  operate  as  designed  if  the  pressure 
within  the  reaction  space  is  high  enough  to  cause  oxidizer 
to  bubble  under  the  lower  edge  of  the  wick.  This  can  occur 
at  high  heat  flux,  where  the  temperature  difference  across 
the  wick,  and  thus  the  pressure  difference,  is  large.  It 
is  likely  that  the  bubbles  of  gas  would  be  reacted  as  they 
flow  to  the  surface  of  the  liquid  [9],  however,  no 
experience  has  been  obtained  with  operation  at  these 
conditions,  and  for  the  present  this  must  be  considered 
to  be  a  failure  mode.  This  limit  is  referred  to  as  the 
liquid  displacement  limit. 

3.  The  system  will  fail  to  operate  as  designed  if  there  is  an 
excessive  temperature  drop  between  the  wick  and  the  load. 
This  results  in  excessive  wick  temperatures,  potentially 
leading  to  failure.  A  variety  of  heat  transfer  limits  are 
encountered,  depending  upon  the  design  of  the  system, 
including:  the  thermal  resistance  of  the  wick,  the  vapor 

transport  resistances,  and  the  condensation  resistance. 

The  condensation  process  is  also  affected  by  the  presence 
of  noncondensibles,  providing  an  allowable  contamination 
limit . 

A  theoretical  model  which  treats  these  aspects  of  the  operation 
of  a  reactive-heat-pipe  combustor  is  presented  in  the  following. 

4. 1  Wick  Analysis 

Figure  13  is  a  sketch  of  the  wick  configuration  to  be  modelled, 
indicating  the  notation  of  the  analysis.  The  wick  has  a  length  (H-L) 


above  the  liquid  surface.  A  liquid  surface  displacement  L  exists  to 
compensate  for  the  pressure  difference  across  the  wick. 

In  cases  of  laminar  combustion  with  free  convection,  the  rate  of 
reaction  is  approximately  proportional  to  the  fourth  root  of  wick 
length  [10].  This  is  not  a  strong  variation,  and  since  the  flow 
pattern  within  the  reaction  space  of  the  combustor  is  not  known,  we 
assume  that  the  burning  rate  of  lithium  is  uniform  over  the  wick 
surface.  The  rate  of  evaporation  of  lithium  from  the  outside  surface 
of  the  wick  is  also  taken  to  be  uniform.  The  system  is  assumed  to  be 
in  a  steady  state  condition  and  all  the  heat  generated  by  reaction  is 
assumed  to  be  transferred  through  the  wick,  in  the  region  above  the 
liquid  surface.  The  wick  heat  flux  q"  is  prescribed,  since  this 
parameter  is  a  design  variable. 

These  assumptions  are  conservative,  since  their  effect  is  to 
increase  the  liquid  pumping  requirements  of  the  wick.  The  circulation 
pattern  within  the  reaction  space  is  expected  to  cause  the  flow  to 
pass  upwards  along  the  wick,  so  that  the  rate  of  reaction  will  be 
greater  towards  the  bottom  of  the  wick,  reducing  the  height  to  which 
the  reacting  fuel  must  be  pumped.  Reaction  energy  produced  toward 
the  bottom  of  the  wick,  in  the  region  H-L  <_  y  H,  is  transferred 
directly  to  the  bath,  causing  lithium  to  evaporate  from  the  liquid 
surface.  Since  the  analysis  neglects  this  energy  removal  mechanism, 
the  lithium  flow  for  evaporation,  and  the  temperature  gradients  in  the 
wick,  are  somewhat  overestimated. 

The  total  mass  flux  of  fuel  for  both  combustion  and  evaporation, 
per  unit  exposed  area  of  wick  is 


Since  the  inside  surface  of  the  wick,  of  length  L,  also  contributes 
to  the  reaction  energy,  these  mass  fluxes  can  be  expressed  in  terms 
of  the  wick  heat  flux  and  the  thermodynamic  properties  of  lithium, 
as  foil ows : 


The  total  lithium  consumption  per  unit  length  along  the  lower 
circumference  of  the  exposed  portion  of  the  wick  is 


*i-L  *  <a-L>  <”c  +  4E}  « 

At  typical  combustor  temperatures,  the  reaction  energy  of 

lithium  and  SF  is  12.5  kW-hr/kg  of  lithium  (4.75  kW-hr/kg  of  SF, 

from  Reference  [3]).  Since  h^  does  not  vary  very  much  in  the 

present  temperature  range,  1100-1300  K,  the  ratio  of  Q  /h  can  be 

r  fg 

taken  as  a  constant.  At  1200  K,  h  =  20818  KJ/K  [11],  yielding 

o  8 

a  -  Qr/hfg  =  2.195. 

Equation  (4.3)  can  then  be  written 

™H_L  "  P  I  (1-3)  (l-Ha-3)  ]  (4. 

xr 

where  8  =  L/H. 

The  lithium  flow  in  the  wick  varies  in  the  vertical  direction 
as  follows 


m'  *  ym^_L/(H-L)  ;  H-L  y  >_  0 


(4. 


The  hydrostatic  and  viscous  pressure  drops  within  the  wick  are 


provided  by  the  capillary  pumping  force 


Aps  -  ApB  +  Apv 


To  calculate  the  total  viscous  pressure  loss,  Darcy's  Law  is 
applied  [12,  15,  16] 


dPy  ^  Ufi’PDw 
dy  "  PfKA„ 

The  total  viscous  pressure  drop  is  given  by 


H-L 
■  / 


Substituting  Equation  (4.5)  into  Equation  (4.7)  and  completing  the 
integration  yields 

U.H2(l-6)2q" 

4pv  *  2pfKtqr  '  <ct+1‘6)  (i 

where 


Assuming  constant  liquid  density,  the  hydrostatic  pressure  drop 


ApB  =  pfg(H-L) 


The  capillary  pumping  force  of  the  wick  is  [12,  15,  16] 


RC  ra 


Since  the  meniscus  radius  at  position  B  is  always  infinitely 

2<r 

large,  we  get  p,gL  =  •£—  and  we  can  write  Equation  (4.11)  as 

f  ra 

follows 

2q" 

APS  "  R-  -  Pf3L  (*.: 

c 

The  capillary  pumping  force,  Apc,  can  have  its  maximum  value  when  R 
reaches  its  minimum  value.  This  occurs  when  [14] 

Rc  "  Rm  “  ($+d')/2  (4 

At  the  maximum  pumping  height  condition,  the  maximum  capillary 
force  just  balances  the  hydrostatic  and  viscous  pressure  loss, 
yielding 


^  -  pfgH  +  q"xH2 

in 


where 


(1-6) 

X  -  -  (a+1-6) 


2pfKtQr 


Solving  for  the  maximum  pumping  height,  we  obtain 


-Pfg  + 


H 
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If  H  is  specified,  the  maximum  heat  flux  is  given  by  Equation 

) 


2a/Rm-prgh 


The  permeability  of  the  wick  can  be  calculated  from  [17] 


K  =  d,2e,3/[122(l-e')2]  ( 

where 

£ 1  -  l-(iTFNd74)  ( 

The  crimping  factor,  F,  has  a  value  of  1.05  [17]. 

In  order  to  calculate  H  ,  we  have  to  know  the  value  of  L  in 

max 

8.  L  can  be  determined  by  noting  that  the  liquid  head  must  equal 
the  pressure  difference  across  the  wick. 

Pwi  -  ”wo  ‘  SPfL  ( 

The  temperature  drop  through  the  wick  is 

where  k  Is  the  effective  thermal  conductivity  of  the  wick.  This 
w 

parameter  is  determined  as  follows  [12] 

\  =  (B’-e)kf/(8’+€)  < 

where 

8'  -  (kf  +  kg)/(kf  -  kg) 

and  £  is  the  volume  fraction  of  the  solid  phase. 

The  temperature  jump  due  to  evaporation  results  from 
nonequilibrium  at  the  liquid  surface,  and  can  be  represented  as 
follows  [13] 
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(4.22) 


where  a 


2RT  /M 
g 


U  ? >v  represents  the  ratio  of  the  bulk  vapor 

/w  \  ^ 


velocity  to  a  characteristic  molecular  velocity.  Equation  (4.22) 
results  from  the  assumption  that  a  is  very  small,  due  to  large  heat 
of  evaporation  of  lithium;  and  also  neglects  the  effect  of  the  mean 
velocity  of  liquid  molecules  at  the  phase  interface  [19],  Since 
c?e  *  cr^  *  1.0  [20],  Equation  (4.22)  becomes 


c\l/2 


2M ' 
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(4.23) 
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The  saturated  vapor  pressure  of  lithium  can  be  represented  as 


[18] 


p  -  exp [A- (B/T) ]  (4.24) 

Knowing  the  temperature  of  the  inner  surface  of  the  wick  and  the 
temperature  after  the  evaporation  jump,  provides  the  following 
expression  for  the  liquid  displacement 

L  -  (exp[A-B(l/Twi41/Tg)](exp[B/Tg]-exp[B/Twi])}  (4.25) 

As  long  as  we  know  L,  by  trial  and  error  we  can  calculate  H 

*  max 

from  Equation  (4.15).  If  H  is  fixed,  we  can  calculate  the  q"  bv 

max 

trial  and  error  from  Equations  (4.16),  (4.20),  (4.23),  and  (4.25). 

When  the  heat  flux  q"  is  small,  $  can  be  set  equal  to  zero  without 
a  large  error,  yielding  a  straightforward  calculation.  The 


thermophysical  properties  required  in  the  calculation  are  discussed 
in  Section  4.4. 


4.2  Transport  Tube  Analysis 

In  order  to  Simplify  the  analysis,  the  flow  of  vapor  to  the 
condenser  and  the  condensation  process  itself  were  considered 
independently,  although  in  many  configurations  flow  and 
condensation  occur  at  the  same  time,e.g.,  in  the  baffled  condensed 
used  in  Tests  4-6.  The  model  under  this  assumption  is  illustrated 
In  Figure  14.  Flow  process  is  taken  to  be  adiabatic  with  friction, 
with  the  vapor  entering  the  passage  without  any  loss.  The  gas 
adjacent  to  the  condenser  surface  is  taken  to  be  "well  mixed"  so 
that  the  concentration  of  noncondensible  gas  is  uniform  over  the 
surface  of  the  condenser. 

The  transport  system  is  composed  of  a  converging  nozzle  and  a 

tube  with  constant  friction  factor  and  constant  cross-sectional 

area.  Adiabatic  vapor  flow  in  the  tube  is  assumed.  Since  the 

bulk  vapor  velocity  leaving  the  wick  is  relatively  small,  e.g., 

2 

0.8  m/sec  at  1000  kW/m  wick  heat  flux  and  1200  K,  we  can  reasonably 
assume  that  the  vapor  In  the  chamber  is  at  rest  without  much  error. 
We  also  assume  that  the  acceleration  process  in  the  inlet  nozzle  is 
isentropic,  and  any  pressure  recovery  as  the  flow  enters  the 
condenser  area  is  neglected.  The  required  pressure  ratio  of  the 
inlet  and  outlet  of  the  nozzle  is  given  implicitly  as  follows  [21] 


CONDENSATION  JUMP 


CONDENSATE  RETURN 


The  Mach  number  and  the  temperature  at  the  inlet  of  the  constant 


area  section  is  then  given  by 

I  2  ,1  n  r/p1i1"Y/Y  1 

k)  -l] 


T±  -  tq/[1  +  (Y-DM^/2] 


The  Mach  number  at  the  outlet  of  the  transport  tube  is  determined 
from  [21] 


and  the  pressure  ratio  and  temperature  ratio  across  the  tube  are 


T  1+(Y-1)M.  /2 

e  i 

-  53  i  ■  ■  ■  i 

Ti  1+(Y-1)M  2/2 
e 


Pe  M.  1+(y-1)M1‘/2 


Me  [l+(Y-DMe2/2 


Given  the  inlet  conditions  to  the  transport  tube,  and  the 
heat  flow  per  unit  area,  Equation  (4.26)  can  be  solved  to  yield 
Pi//po*  Substituting  the  value  of  P^/po  into  Equations  (4.27)  and 
(4.28),  then  provides  the  Mach  number  and  temperature  at  state  i. 

The  quantity  fLt/D  is  taken  as  a  parameter  in  the  analysis,  knowing 
this  quantity  and  the  conditions  at  state  i,  simultaneous  solution 
of  Equations  (4.29)  and  (4.30)  yields  the  Mach  number  and  temperature 
ratio  at  the  exit  of  the  tube.  Equation  (4.31)  finally  provides  the 
condenser  pressure. 


As  long  as  the  flow  is  not  sonic  at  the  exit  of  the  transport 
tube,  solution  of  Equations  (4 . 26 >-(4 . 31)  yields  the  pressure  in  the 
region  adjacent  to  the  condenser  [21].  If  the  flow  is  choked,  the 


pressure  in  the  condenser  region  is  determined  by  the  cooling 
characteristics  of  the  condenser  and  there  is  a  pressure 
discontinuity  at  the  exit  of  the  transporting  tube.  The  choked 
condition  is  defined  by  M  ■  1,  and  the  corresponding  inlet  Mach 
number  is  given  implicitly  by 


(Y+l)Mi2/2 

l+(Y-l)Mi2/2 


(4.3 


In  order  to  return  the  condensate  in  the  condenser  space  to 
the  combustor  chamber,  a  certain  pressure  head  must  be  applied  to 
the  condensate  to  compensate  for  the  pressure  difference  between  the 
combustion  chamber,  and  the  viscous  loss  in  the  return  tube.  The 
loss  in  the  return  tube  can  be  calculated  by  conventional  methods. 


4.3  Condenser  Analysis 

4.3.1  Condensation  in  the  Absence  of  Noncondensibles.  If  there 
is  no  noncondensible  gas  in  the  condenser,  the  thermal  resistances 
from  the  vapor  space  to  the  inner  surface  of  the  condenser  are  the 
resistance  of  the  condensate  film  on  the  heat  exchanger  surface, 
and  the  interfacial  resistance  between  the  vapor  and  liquid  film. 

For  the  liquid  film,  we  assume  that  subcooling  of  the 


condensate  may  be  neglected  due  to  the  high  thermal  conductivity  of 
liquid  lithium.  The  temperature  variation  is  illustrated 
schematically  in  Figure  15.  The  condensation  resistance  due  to  the 
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liquid  film  is  [27] 


T..-T  . 
f  i  wi 


3  1] 

BPf(prpg)kf  ht« 


(4.33) 


In  Equation  (4.33),  a  correction  factor  of  1.2  [27]  is  used  to  modify 

the  theoretical  value  for  the  effect  of  the  downward  flowing  vapor  in 

reducing  the  thickness  of  the  condensate  film  and  the  mixing  action 

of  ripples  in  this  film.  This  correlation  is  based  on  the 

evaluation  of  the  properties  at  the  mean  film  temperature 

(T-.+T  .)/2,  while  the  latent  heat  h.  is  to  be  taken  at  T_.  [28], 
n  wi  rg  ii 

If  the  film  is  turbulent,  then 


T_.-T  . 
f  i  wi 


-  0.0134 


3  2  \1/3/0  2  3  ,0,4 

kf  Pf  8  Pf  gLc 


(4.34) 


For  the  cases  considered,  the  film  is  always  in  a  laminar  region. 

2 

Considering  the  extreme  case,  qj]  =  1000  KW/m  and  Twi  =  1400  K,  the 

film  is  laminar  when  L  is  below  0.6m. 

c 

The  interfacial  resistance  is  essentially  the  reciprocal  of  the 
evaporative  resistance  [13] 


*c  r 2m1  /  rv  pfi 

hfg  -  (_  ttR J  -  ^172 


(4.35) 


where  pressures  are  related  to  the  temperatures  through 


p  =  exp[A-(B/T) ] 


4.3.2  Condensation  in  the  Presence  of  Noncondensibles.  When  a 


noncondensible  gas  is  present,  there  is  a  gas  film  adjacent  to  the 
liquid  film.  The  heat  transfer  across  the  gas  film  is  the  sum  of  the 
sensible  heat  transferred  across  the  layer  and  the  heat  of 
vaporization  of  lithium  vapor  that  diffuses  across  the  layer  and 
condenses  on  the  liquid  film.  The  basic  equation  for  the  gas 
film  is 


q"  =  h(T  -T  , )  +  K  Mh  (p  -p  ,) 
g  g'  g  fg  g  8 


(4.3l 


If  no  forced  convection  is  present,  the  convective  effect  is  measured 
by  the  Grashof  number 


Gr  *  8Lc3(‘V-ps>/V''2 

The  mass  transfer  coefficient  is  given  by  [22] 


(4.3 


K  L  RT  p 
g  c  am 


1.02(Gr  •  Sc) 


0.373 


(4.3 


Li-Ar 
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for  10  <  Gr  <  10  .  In  this  equation,  p^  is  the  logarithmic  mean 


of  the  partial  pressures  of  noncondensible  gas  at  the  interface  and 
in  the  bulk  mixture. 


=  <Pg-Pg.)/M(P-Pg.>/(P-Pg)] 


(4.31 


The  sensible  heat  transfer  coefficient  is  obtained  from  [28] 


hL 
_ c 

K 


0  25 

0.55  (Gr  •  Pr) 


(4.4i 


for  10 


g 

4 


<  Gr  Pr  <  10 
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Since  the  rate  of  mass  transfer  is  less  than 


55 


2  2 

0.05  kg/m  s,  when  the  heat  flux  in  the  condenser  is  1000  kW/m  ,  we 

can  simply  use  the  coefficient  given  in  Equation  (4.40)  without 

introducing  much  error  due  to  the  bulk  mass  transport  effect  [13]. 

Since  interfacial  conditions  are  not  known,  Equations  (4.36)- 
(4.40)  must  be  solved  by  trial  and  error.  The  properties  of  gas 
mixtures  used  in  the  calculation  are  summarized  in  Appendix  A. 

Since  the  pressure  recovery  of  the  flow  leaving  the  transport 
tube  has  been  neglected,  the  vapor  in  the  gas  space  adjacent  to  the 
condenser  is  somewhat  superheated.  Since  the  enthalpy  of  superheating 
is  small  compared  to  h^,  however,  this  contribution  has  been 
neglected.  Superheating  improves  heat  transfer  characteristics 
[23],  therefore,  neglecting  this  phenomenon  makes  the  present 
calculations  somewhat  conservative. 


4.4  Thermophysical  Properties 

The  properties  of  lithium  and  argon  are  mainly  selected  from 
References  [11,  18].  Although  the  reported  properties  are  not  quite 
the  same  for  these  two  sources,  the  discrepancies  are  not  large 
over  the  present  range  of  calculations.  The  atomic  weights, 
melting  points,  boiling  points,  and  critical  points  of  lithium  and 
argon  are  summarized  in  Table  3. 

4.4.1  Properties  of  Lithium.  The  detailed  properties  of  lithium 
are  presented  in  Table  4  [11].  For  convenience  in  the  calculations 
the  properties  of  liquid  lithium  can  be  correlated  by  the  following 
relation  [18] 


Atomic  Weight  (kg/kg-mole) 

6.940 

39.944 

Melting  Point  (*K) 

453.7 

83.78 

Boiling  Point 

at  760  mm  Hg  (*K) 

1615 

87.29 

Critical  Temperature  (*K) 

3800 

150.86 

Critical  Pressure  (bar) 

970 

50 

Density  at  Critical 
Point  (kg/m^) 


100 


536 


Table 


1.162x10 


4  (Continued) 


interpolation  in  Reference  [11] 


a  =■  0.447  -  1.07  x  10  T  -  1.351  x  10  T  (4.42 
k  =  21.874  +  0.056255  T  -  1.9325  x  10~5  T2  (4.43 
p  -  exp  [23.06  -  18569. 20/T]  (4.44 
M  -  exp  [-7.09  +  132.70/T  -  1.190  x  10-3  T]  (4.45 


3  2 

where  T  is  in  K,  p  in  kg/m  ,  0  in  N/m,  k  in  W/mK,  p  in  N/m 

2 

and  p  in  Ns/m  . 

Since  similar  correlations  for  lithium  vapor  are  not  available, 
the  data  in  Reference  [11]  was  correlated  to  obtain  the  following 


formulas  for  saturated  liquid  lithium  vapor. 

Ml±  *  exp  [2.02  -  1.4  x  10~4  (1400-T) ]  (4.46 

log1()  p  «  -  49.052  +  14.988  logl0C-589  +  2.0275  T 

-  4.375  x  10"4  T2)  (4.47 

k  -  5.641  x  10~2  T1/2  -  6.28  x  10_4  T  -  1.05975  (4.48 

V  -  1.29  x  10'5  +  10'8  x  (T-1300)  (4.49 


The  temperature  ranges  for  the  above  correlations  are  summarized 
in  Table  5. 

4.4.2  Properties  of  Argon.  From  the  data  of  Reference  [11], 
it  is  found  that  most  of  the  properties  of  argon,  except  density,  are 
not  sensitive  to  pressure  changes.  For  example,  the  viscosity  at 
900  K  changes  only  about  one  percent  when  the  pressure  is  changed 


Table 


from  1  to  80  bar.  Therefore,  the  properties  of  argon  are  approximated 
by  the  values  at  a  pressure  of  1  bar,  and  are  presented  in  Table  6 
[11].  The  thermal  conductivity  of  argon  in  the  range  1000-1400  K, 
at  1  bar,  can  be  correlated  by  the  following  equation 

k  -  4.76  x  10~2  +  2.6  x  IQ-5  (T  -  U50)  (4.50) 

and  the  viscosity  in  the  range  of  1000-1500  K  at  1  bar  can  be 
correlated  by 

y  -  6.679  x  10“5  +  3.5  x  10-8  (T  -  1300) 


(4.51) 


Temperature 


K 


Thermoconductivity 
(x  103) 

W/m-K 


Viscosity 
(x  107) 
Ns/m3 


Specific 

Heat 

(Cp) 

KJ/kg-K 


CHAPTER  V 


RESULTS 

The  analytical  results  for  the  wick,  and  for  the  vapor 
transporting  and  condensation  process,  will  be  discussed  separately. 
This  will  be  followed  by  consideration  of  the  complete  system  for 
some  typical  conditions. 

5.1  Wick 

Two  thermal  resistances  are  considered  here,  the  wick  itself  and 
the  evaporative  interface.  Figure  16  is  an  illustration  of  the  heat 
transfer  characteristics  of  a  wick  consisting  of  two  layers  of  100- 
mesh  screen,  which  is  the  arrangement  employed  in  the  tests.  The 
total  temperature  drop  increases  almost  linearly  with  increasing 
heat  flux.  Since  the  thermal  conductivity  of  liquid  lithium  varies 
slightly  with  temperature,  the  resistance  of  the  wick  is  relatively 
independent  of  temperature.  The  evaporative  resistance  increases  as 
the  mean  temperature  level  decreases. 

As  shown  in  Figures  17  and  18,  increasing  the  number  of  layers 
of  screen,  and  reducing  the  mesh  size,  increases  the  thermal 
resistance  of  the  wick,  however,  the  evaporative  resistance  only 
depends  on  the  mean  temperature  level  and  the  heat  flux.  The  linear 
relationship  between  the  temperature  difference  and  the  wick  heat 
flux  is  also  valid  in  Figures  17  and  18.  The  wider  range  of  the 
temperature  difference  can  be  conveniently  plotted  with  a  semi- 
logarithmic  scale. 


Figure  18  Wick  Thermal  Resistance  with  Different  Mesh  Sizes 


Increasing  the  mesh  will  increase  the  capillary  ability  of  the 

wick.,  but  it  will  also  increase  the  viscous  loss.  Therefore,  there 

is  an  optimum  mesh  for  any  given  operating  condition.  Figure  19 

is  an  illustration  of  the  pumping  height,  H,  and  the  net  pumping 

height  H-L,  for  various  mesh,  sizes  for  heat  fluxes  of  100  and 

2 

1000  kW/m  ,  at  1200  K.  We  can  see  that  the  maximum  pumping  height 
occurs  for  approximately  a  75  mesh,  wick,  under  the  specified 
conditions . 

Figures  20  and  21  are  illustrations  of  the  variation  of  pumping 
height  and  net  pumping  height  as  a  function  of  heat  flux,  with 
temperature  as  a  parameter  for  mesh  size.  The  net  pumping  height 
always  decreases  with  increasing  heat  flux,  and  is  generally  lower 
at  higher  operating  temperatures.  The  pumping  height  itself  reaches 
a  minimum  with  increasing  heat  flux,  and  is  generally  higher  at 
higher  temperatures .  The  minimum  occurs  since  L  becomes  larger  at 
higher  heat  fluxes,  reducing  the  height  that  the  lithium  must  be 
pumped  through  the  exposed  portion  of  the  wick,  which  compensates 
for  the  increased  lithium  flow  rate.  The  temperature  effect  results 
from  the  reduced  density  of  lithium  at  higher  temperatures,  which 
also  increases  L  and  reduces  pumping  requirements  in  the  exposed 
portion  of  the  wick. 

The  net  pumping  heights  illustrated  in  Figures  19-21, are 
relatively  small,  which  limits  the  flexibility  of  design  since 
reasonable  wick  heights  can  be  necessary  to  provide  sufficient 
reaction  area  within  a  given  volume.  Larger  net  pumping  heights 
can  be  obtained  by  employing  an  arterial  wick.  In  this  case,  the 


WICK  HEAT  FLUX  (KW/M1*) 


Figure  20  Wick  Pumping  Heights  for  Different  Temperatures 


liquid  flows  through  an  independent  tube,  and  is  not  just  confined 
to  the  cross-section  of  the  wick.  Only  the  final  distribution  of  the 
liquid  across  the  region  between  arteries  occurs  within  the  wick 
itself.  In  this  case,  the  viscous  pressure  drop  in  the  wick  can  be 
made  much  smaller  than  in  the  conventional  wick  arrangement,  which 
increases  the  net  pumping  height. 

Table  7  summarizes  the  pumping  heights  for  four  different  mesh 
sizes  using  an  arterial  wick.  When  an  arterial  wick  is  used,  finer 
mesh  screen  is  favored,  and  the  pumping  height  is  independent  of  the 
heat  flux.  The  pumping  heights  are  also  substantially  greater 
than  the  results  shown  in  Figure  19,  where  the  liquid  flow  is 
confined  to  the  wick  cross-section. 


5.2  Vapor  Transport  and  Condensation 


The  losses  in  this  section  include  the  vapor  flow  resistance 
and  resistances  in  the  condenser.  The  condenser  resistance  involves 
the  resistance  of  the  film,  the  resistance  at  the  interface,  and  the 
resistance  of  diffusion  through  any  noncondensible  gas  layer.  In 
order  to  highlight  losses  in  the  vapor  transport  tube,  the  area  of 
the  condenser  is  taken  to  be  ten  times  larger  than  that  of  the  wick, 
while  the  vapor  transport  tube  is  taken  to  be  only  .1  that  of  the 
wick.  In  order  to  fix  ideas,  the  combustor  geometry  used  in  the  tests 
employed  A  /A  =0.23;  and  fL/D  in  the  range  0-10,  depending  upon 
operation  conditions  and  the  appropriate  selection  of  transport  tube 
length  (in  view  of  the  fact  that  vapor  flow  and  condensation  occur 
simultaneously  in  a  portion  of  the  arrangement  the  length  can  only 


be  approximated) .  Figure  22  is  an  illustration  of  a  portion  of  the 


Table  7 


for  Arterial  Wicks  (mm) 


Temperature  (K) 


Mesh 

Size 

1100 

1250 

30  x 

30 

334 

322 

309 

50  x 

50 

556 

537 

515 

100  x 

100 

1113 

1074 

1031 

250  x 

250 

2782 

2685 

2577 

74 


results  obtained  for  the  complete  vapor  transport  system  and 
condenser  in  the  absence  of  noncondensibles.  Figure  23  shows 
similar  results  when  a  noncondensible  gas  is  present.  The 
reference  temperature  shown  on  these  figures  is  the  temperature 
at  the  inlet  of  the  vapor  transport  tube. 

In  general,  the  temperature  drop  for  all  steps  in  the  heat 
transport  system  increases  as  the  heat  flux  increases,  and  as  the 
mean  temperature  of  the  vapor  decreases.  The  latter  effect  is  the 
result  of  the  lower  vapor  pressure  and  lower  thermal  conductivity  of 
liquid  lithium  at  these  conditions. 

In  Figure  22,  the  results  for  various  values  of  fL^/D  are 
identical  at  low  heat  fluxes,  however,  the  temperature  drop 
increases  more  rapidly  for  higher  values  of  fLt/D  at  the  higher  heat 
fluxes.  The  results  for  fLt/D  =  10  at  1100  K  illustrate  the  phenomena 
of  choking  in  the  vapor  transport  tube.  As  the  heat  flux  increases, 
a  condition  is  reached  where  the  temperature  drop  increases  very 
rapidly  with  increasing  heat  flux.  Eventually  a  condition  is 
reached  where  the  heat  flux  reaches  a  maximum  value,  irrespective 
of  the  available  temperature  difference,  for  the  given  inlet 
temperature.  At  this  condition,  the  vapor  transport  tube  can  only 
deliver  the  maximum  heat  flux,  irrespective  of  conditions  in  the 
condenser,  and  the  flow  at  the  exit  of  the  transport  tube  is 
choked.  The  maximum  heat  flux  under  choked  conditions  decreases 
as  the  friction  parameter  is  increased. 

Operation  in  the  choked  region  may  occur  during  startup  of  the 
system  when  the  condenser  is  at  low  temperatures.  Through  choking, 


Figure  23 


Flow  Losses  in  the  Vapor  Transport  System  with  Argon 
in  the  Condenser 


only  a  limited  amount  of  heat  can  be  supplied  to  the  load,  while 
the  remainder  of  the  energy  of  reaction  serves  to  increase  the 
temperature  level  of  the  wick  area.  As  the  temperature  of  the 
wick  increases,  the  choked  heat  flux  increases  as  well,  as  shown 
in  Figure  22 .  Therefore ,  more  energy  is  transported  to  the  load 
as  the  design  wick  temperature  level  is  approached,  and  choking  ceases 
to  be  a  limiting  factor.  For  the  arrangements  used  during  the  test 
program  discussed  earlier,  fLfc/D  is  small  and  losses  in  the  vapor 
transport  tube  were  negligible,  except  during  startup. 

In  Figure  23,  the  reference  temperature  has  been  transferred 
to  the  condenser  area,  and  noncondensibles  are  considered.  There 
is  no  choking  in  the  vapor  transport  system  for  the  temperature 
levels  shown  on  the  figure.  The  thermal  resistance  due  to  the 
noncondensible  gas  is  sufficiently  large  so  that  the  total  heat  flux 
in  the  system  is  not  sufficient  to  cause  choking.  It  is  evident 
that  even  a  low  level  of  noncondensibles  can  provide  a  thermal 
resistance  equivalent  to  that  of  the  film  condensation  process 
itself.  Higher  levels  of  noncondensibles  will  be  considered  in  the 
next  section. 

5.3  Thermal  Resistances  in  Some  Typical  Systems 

Temperature  distributions  within  some  typical  system  are 
illustrated  in  Figures  24  and  25.  Figure  24  considers  operation 
when  there  is  no  noncondensible  gas  in  the  system,  while  Figure  25 
considers  the  effect  of  noncondensibles.  Both  figures  are 
constructed  for  a  condenser  wall  temperature  of  1200  K  and  a  heat 


flux  of  1000  kW/m^. 


Figure  24  Temperature  Distributions  in  a  Reactive-Heat-Pipe 
with  No  Argon  in  the  Condenser 


The  baseline  system  in  Figure  24  is  represented  by  curve  1, 
which  is  constructed  for  a  case  where  the  wick  area  equals  the 
condenser  area,  the  area  of  the  transport  tube  is  .1  that  of  the 
wick,  and  the  frictional  resistance  of  the  transport  tube  is 
negligible.  In  this  case,  the  total  temperature  drop  is  31.5  C. 

The  major  resistances  in  the  system  are  the  temperature  drop  through 
the  wick  and  the  temperature  drop  caused  by  acceleration  of  the 
flow  into  the  transport  tube. 

Curve  2  in  Figure  24,  involves  increasing  the  condenser  area 
by  a  factor  of  10.  This  results  in  only  a  small  reduction  in  the 
total  temperature  drop,  since  the  condenser  resistance  is  not 
a  controlling  factor  for  the  baseline  system.  Curve  3  involves  a 
finite  resistance  in  the  vapor  transport  tube;  the  effect  here  is  to 
increase  the  temperature  drop  across  the  transport  tube,  resulting 
in  an  overall  temperature  drop  to  the  system  of  38  K.  Therefore,  it 
is  evident  that  the  design  of  various  aspects  of  the  system  can 
shift  the  controlling  heat  transfer  resistance  from  wick  to  vapor 
transport  tube  to  condenser. 

Figure  25  compares  the  same  baseline  system,  with  conditions 
which  have  various  levels  of  noncondensibles  in  the  condenser. 

Curve  1  represents  the  baseline  system,  with  its  temperature  drop 
of  31.5  C.  When  argon  is  present  in  the  system  with  a  mole  fraction 
of  .1  in  the  condenser  area  (curve  5),  the  overall  temperature  drop 
in  the  system  increases  to  372  C.  Reducing  the  argon  mole  fraction 
to  .001  (curve  3),  only  reduces  the  temperature  drop  to  262  C. 
Therefore,  the  effect  of  noncondensibles  is  very  nonlinear  and  a 


substantial  resistance  is  provided  by  very  small  quantities  of 
noncondensibles.  Of  all  the  loss  mechanisms  examined  during  this 
study,  it  appears  that  the  presence  of  noncondensible  gases  in  the 
condenser  is  one  of  the  most  important  factors. 

There  are  several  ways  of  overcoming  the  noncondensible  problem. 
One  approach  involves  reducing  the  heat  flux  on  the  condenser  surfaces 
For  example,  the  results  in  Figure  25  for  curve  2,  indicate  that  the 
total  temperature  drop  is  reduced  to  57  C,  for  a  mole  fraction  of 
.001,  when  the  condenser  area  is  increased  by  a  factor  of  10. 

For  many  types  of  loads,  particularly  in  the  case  of  gas  cycles, 
limitations  on  the  allowable  heat  exchanger  flux  will  result  in 
conditions  of  this  type.  Another  approach  involves  the  use  of 
baffles  or  constricted  vapor  flow  passages,  to  increase  the  vapor 
velocity  past  the  condenser  surface.  This  increases  the  surface 
heat  and  mass  transfer  coefficients  in  the  noncondensible  layer 
above  the  natural  convection  calculation  values  used  in  the  present 
calculations.  Baffles  also  tend  to  sweep  the  noncondensible  material 
away  from  the  surface  of  the  condenser.  Both  of  these  effects 
reduce  the  heat  transfer  resistance  due  to  noncondensihles .  When 
all  these  techniques  are  used,  in  conjunction  with  practical  levels 
of  control  of  noncondensibles  in  the  system,  acceptable  designs 


should  be  feasible. 


CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 


6.1  Summary 

The  overall  objective  of  the  present  study  was  to  investigate, 
both  analytically  and  experimentally,  the  characteristics  of  a 
modified  liquid  metal  combustor  employing  the  lithium-sulfur 
hexafluoride  reactant  combination.  The  specific  objectives  of  this 
study  were  as  follows: 

1.  Design  and  test  a  small  scale  version  of  the  reactive 
heat  pipe  combustor  in  order  to  observe  its  operation. 

2.  Complete  a  theoretical  analysis  of  the  reactive  heat  pipe 
syst '  Ln  order  to  obtain  the  design  characteristics  of 
the  wick,  and  the  heat  transport  system. 

The  oxidizer  was  injected  into  the  space  inside  the  free-standing 
wick,  and  reacted  with  lithium  vapor  evaporated  from  the  wick. 

Lithium  vapor  evaporating  from  the  outer  surface  of  the  wick  served 
to  transfer  thermal  energy  to  the  load.  The  cylindrical  combustor 
wall  above  the  liquid  level  acted  as  the  condenser  surface  of  the  heat 
pipe.  Heat  loss  from  the  combustor  through  the  surface  was  controlled 
by  sweeping  them  toward  a  collector,  through  a  baffle  system  located 
parallel  to  the  condenser  surface. 

The  system  was  operated  in  the  wall  temperature  range  1088-1240 

2 

K  with  wick  power  densities  in  the  range  0-260  kW/m“.  No  corrosion 
problems  were  observed  using  type  316  stainless  steel  for  the 
components  of  the  combustor.  The  capability  of  the  system  to  operate 


under  variable  load,  including  idling,  as  well  as  restarting,  was 
also  demonstrated. 

The  oxidizer  flow  rate  was  varied  in  order  to  accommodate 
different  load  conditions.  At  a  fixed  wall  temperature,  the 
temperature  and  pressure  in  the  reaction  space  did  not  vary 
significantly  as  the  load  was  changed.  The  condenser  wall 
temperatures  were  nearly  constant  at  full- load  conditions,  which 
demonstrates  the  desirable  characteristics  of  the  heat  pipe 
configuration.  A  major  practical  problem  encountered  during  the 
testing  was  the  presence  of  small  quantities  of  noncondensibles  in 
the  system.  The  noncondensibles  acted  to  increase  the  heat  transfer 
resistance  of  the  condenser.  Problems  were  circumvented  by  using 
a  noncondensible  collection  tube  and  installing  baffles  in  the  system. 

The  theoretical  analysis  considers  screens  of  two  and  four  layers 

with  mesh  sizes  ranging  from  30  to  250.  System  temperatures  were 

varied  from  1100-1400  K,  and  the  heat  flux  through  the  wick  was  in 

2 

the  range  0-1000  kW/ra  .  The  presence  of  "oncondensibles ,  as  argon, 
were  considered  in  the  mole  fraction  range  0-.1.  The  vapor  in  the 
condenser,  was  conservatively  considered  to  condense  under  free- 
convection  conditions.  The  pressure  drop  characteristics  in  the 
vapor  transport  system  between  the  wick  and  the  condenser  were  also 
considered. 

6 . 2  Conclusions 

The  major  conclusions  of  the  study  are  r.s  follows: 

1.  For  a  given  heat  flux,  the  maximum  pumping  height  can  be 

increased  by  using  more  layers  of  screen  in  order  to  reduce 
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the  viscous  losses.  There  is  an  optimum  mesh  size  for 
given  conditions  which  yields  the  maximum  pumping  height; 
below  this  size,  reduced  capillary  forces  limit  performance; 
above  this  size,  increased  viscous  forces  limit  performance. 

2.  The  liquid  displacement  at  the  base  of  the  wick  increases 
as  the  thickness  of  screen,  temperature,  and  heat  flux  of 
the  wick  is  increased.  This  effect  can  limit  the  use  of  a 
large  number  of  screens  in  order  to  reduce  viscous  losses. 

3.  The  thermal  resistance  of  the  evaporation  process,  the 
vapor  transport  process  and  the  condensation  process, 
increases  as  the  temperature  level  in  the  system  decreases. 
At  low  temperature  levels  choking  can  occur  in  the  vapor 
transport  system  which  limits  the  energy  delivery  to  the 
load. 

4.  The  most  significant  thermal  resistance  appears  to  be  due 
to  the  presence  of  noncondensibles  in  the  system.  The 
noncondensibles  can  be  controlled  by  limiting  condenser 
heat  fluxes,  using  a  baffle  to  sweep  noncondensibles  from 
the  heat  exchanger  surface  [13],  and  increasing  vapor 
velocities  along  the  surfaces  of  the  condenser  in  order 
to  reduce  the  dlffusional  resistance  of  the  layer  of 
noncondensib les . 

5.  The  experiments  demonstrated  relatively  uniform  temperature 
distributions  on  a  load  heat  exchanger,  easy  repair  and 
maintenance  of  the  system,  no  danger  of  local  subcooling 

of  the  wick  during  high  heat  loads,  compactness  of  the 


to  load  changes.  Observations  of  the  combustor  after  the 
test  did  not  Indicate  any  substantial  degree  of  corrosion, 
however,  the  present  test  periods  were  relatively  short, 
extending  only  to  four  hours  of  operation.  Longer 
duration  testing  is  desirable,  in  order  to  fully 
evaluate  the  effects  of  corrosion. 

6.3  Recommendations  for  Further  Study 

In  many  practical  applications  higher  wick  pumping  heights  than 
can  be  accommodated  by  the  design  considered  in  this  study  would  be 
desirable.  The  analysis  indicates  that  the  use  of  arterial  wicks  would 
provide  such  capability,  and  the  feasibility  of  systems  of  this  type 
should  be  examined. 

With  regard  to  the  analysis,  the  present  investigation  has  been 
limited  to  natural  convection  conditions  in  the  condenser.  This  is 
the  most  conservative  estimation  of  the  effect  of  noncondensibles  and 
forced  convection  phenomena  should  be  considered  as  well.  The 
analysis  has  also  been  limited  to  the  case  where  the  quantity  of 
noncondensible  in  the  system  is  sufficiently  large  so  that  the 
increased  concentration  at  the  surface  of  the  condenser  does  not 
significantly  deplete  the  ambient  concentration  of  noncondensibles. 
Further  analysis  would  be  desirable  at  the  limit  of  low  levels  of 
non condensibles,  following  the  lines  indicated  by  the  work  of 
Rohsenow  and  Kroger  [29]. 
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APPENDIX  A 


CALCULATION  OF  THE  THERMAL  RESISTANCE  OF  CONDENSATION  FOR  TEST  3 


To  calculate  the  thermal  resistance  of  condensation,  the 
combustor  was  run  at  a  steady  state  condition  so  that  the  heat 
generated  was  equal  to  the  heat  output.  To  simplify  the  calculations 
the  upper  edge  of  the  radiation  shield  Is  taken  to  be  above  the  level 
of  the  liquid  surface,  and  the  top  surface  of  the  combustor  Is 
Insulated.  Because  the  heat  conduction  through  the  Insulation  Is 
small,  we  can  assume  that  all  the  heat  Is  transferred  through  the 
exposed  combustor  wall.  The  thermal  resistance  at  condenser  can 

be  defined  as  follows: 

_  1 


rc  "  h  A 

c  ex 


T  -T  . 
g  cl 


(A.l) 


As  discussed  In  Section  3.2,  we  can  neglect  the  effect  of  heat 
conduction  from  the  condenser  wall  below  the  upper  edge  of  the  shield 
due  to  the  rather  uniform  temperature  distribution.  The  height  of 
the  condenser  wall  above  the  liquid  surface  can  be  directly  measured 
by  vision  or,  more  precisely,  calculated  from  the  properties  of  the 
reactants  and  the  products  as  described  In  Reference  3. 


The  temperature  of  the  Inner  surface  of  the  condenser  is 
I  D 
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where 


Q  -  m  Q 
x  ox  r 


From  Figure  5,  we  can  get  the  temperature  of  the  inner  surface  of  the 
wick  by  measuring  the  combustor  pressure.  Equation  (4.20)  and  (4.23) 
yield  the  temperature  of  the  outer  surface  of  the  wick.  If  the 
radial  vapor  velocity  is  small,  we  can  simply  take  the  temperature 
of  the  outer  surface  of  the  wick  to  be  the  stagnant  temperature 
outside  the  condenser. 

At  80  minutes  from  start-up  in  Test  3,  the  temperature 
distribution  of  the  outer  surface  of  the  condenser  ranged  from  1150 
K  to  1158  K,  the  combustor  pressure  inside  the  wick  cylinder  was  18 
mm  Hg,  the  oxidizer  flow  rate  was  1.32  Kg/hr  and  the  exposed  length 
of  the  condenser  was  115  mm,  which  was  approximately  equal  to  the 
height  above  the  liquid  surface.  The  comparison  of  the  actual  and 
the  expected  thermal  resistances  at  this  instant  of  time  is  computed 
in  the  following. 

From  Figure  5,  we  know  the  inner  surface  temperature  of  the 

wick  is  approximately  1220  K.  The  thermal  conductivity  of  the 

condenser  wall  is  24  W/m-k  [28].  The  wick  heat  flux  is  approximately 
2 

260  Kw/m  .  From  Equation  (A. 2),  the  inner  surface  temperature  of 
the  condenser  is  1189  K.  Following  the  analysis  in  Section  4.1, 
the  temperature  drop  which  includes  the  drop  through  the  wick  and  the 
jump  due  to  nonequilibrium  evaporation  is  3.92  K.  The  bulk  vapor 
velocity  leaving  the  wick  surface  is 


20818x1.494x10 


■  8.36  m/s 

Since  V  is  small,  we  can  take  (1220  -  3.92)K  to  be  the 
temperature  right  before  the  condensation  process.  Then  the 
total  temperature  drop  through  the  condensation  process  is  27  K. 

If  there  was  no  noncondensible  gas  in  the  combustor,  following  the 
analysis  in  Section  4.1,  the  temperature  drop  for  the  condensation 
process  is  1.16  K.  The  resistance  ratio  of  the  experimental  and  the 
expected  value  is  approximately  24. 


4 
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APPENDIX  B 


CALCULATION  OF  PROPERTIES  FOR  GAS  MIXTURES 


B.l  Grashof  Number 

The  Grashof  number  is  defined  as 
L  3(Ap/p)gp2 

Gr  -  — - 1 -  (B. 

U 

In  order  to  calculate  the  density  terms,  the  vapor  gas  mixtures  is 
assumed  to  be  an  imperfect  gas.  The  total  pressure  across  the 
boundary  layer  is  a  constant,  therefore. 


where  Z  is  the  compressibility  factor.  For  condensation  processes, 

which  occur  at  moderately  low  pressures,  we  can  take  Z  , /Z  1 

8  8 

[24],  yielding 


The  molecular  weight  of  the  mixtures  can  be  calculated  as 
follows 


(B. 


M  -  (1  -  W.  )  M.,  +  W.  M.  (B. 

g  Ar,g  Li  Ar,g  Ar 

M  ,  -  (1  -  W.  ,)  M.  .  +  W.  ,  M.  (B. 

g'  Ar ,  g  Li  Ar,g  Ar 

Where  W  and  W.  ,  are  the  mole  fractions  of  noncondensible  gas 
Ar ,  g  Ar ,  g 

in  the  space  far  from  the  liquid  film  surface  and  at  the  surface, 
respectively.  The  mole  fractions  can  be  calculated  from  the 


following  equations 


W.  -  1  - 
Ar,g  p 


■Ll.g' 


Ar,g’ 


The  viscosity  of  argon-lithium  vapor  mixtures  at  low  density 
can  be  well  approximated  by  Wilkes'  semi-empirical  formula  [25]. 

„  WLiULi _  WArUAr  , 

WLi  +  ffAr  *Li-Ar  «Li  V-Li  +  tfAr 
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and  UAr  and  v*Li  are  the  viscosity  of  argon  and  lithium  vapor  at  the 
system  temperature  and  pressure. 

The  densities  of  the  argon-lithium  vapor  mixture  are  determined 


p  .  *  p  H  ,  /RT 

mix  *  mix 


■  ’(“Li  +  MAr  5Ar)/RT 


(B.12: 


Considering  the  changes  of  molecular  weight  of  the  lithium  vapor  due 
to  the  presence  of  diatomic  molecules ,  the  correlation  for  the 
molecular  weight  of  saturated  lithium  vapor  is  applied. 


B.2  Prandtl  Number 


The  Prandtl  number  Is  defined  as 


The  value  of  y  has  already  been  given  in  Equation  (B . 8) .  The 
value  of  is  given  by  a  semi-empirical  formula  similar  to 

Equation  (B.8)  [25]. 


k  WLi  ^i  WAr  kAr 

“lx  "  WLi  +  WAr  ^Li-Ar  WLi  *Ar-Li  +  WAr 


where  <}>. 


Li-Ar 


and  <J>^r_^^  were  given  by  Equations  (B.9)  and  (B 


The  specific  heat  of  argon-lithium  vapor  mixtures  can  t 


calculated  from 


(WT.  C  +  WA  C  )/(WT  .  M.4  +  WA  Ma  ) 
v  Li  p,  .  Ar  pA_  Li  Li  Ar  Ar 


where  C  and  C  are  taken  as  the  average  properties  in  the 
pLi  pAr 

range  1000-1400  K. 


B.3  Schmidt  Number 


The  theoretical  equation  for  the  mutual  defusion  coefficient 
in  a  low  pressure  argon-lithium  vapor  mixture  is  [26] 


Li-Ar 


1.858  x  10"7  T1’5  (l/^Li  +  1/M^)°*5 
p  5L1-Ar  "d 


where  DLi_Ar  is  in  m  /s,  p  in  atm.,  a^_Ar  in  A.  is  the 
collision  integral  and  is  a  function  of  kT/e^.^j.*  vaiues  °i 

this  parameter  can  be  found  in  Table  2.2  of  Reference  [26].  The 


constants  and  aLj_Ar  are  obtained  by  use  of  the  following 


combining  rules 


Li-Ar  =  rLi  "Ar 
'  k  k 


ffLi-Ar  =  2  (aLi  +  aAr) 


The  Lennard-Jones  potential  parameters  a ^  and  e^/k  can  be  found  in 


Table  2.3  of  Reference  [26].  They  are 


a.  =  3.542  A 
Ar 


93.3  K 


The  above  parameters  are  not  available  for  lithium  vapor,  therefore, 
we  use  the  following  approximate  rules  [26] 


—  =»  0.75  (T  )  (B. 

k  c  Li 
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where  T  is  in  K,  V  is  in  cm  /g  mole  and  Or  .  is  in  A.  From  Table 
c  c 

3  we  obtain 


=  2850  K 


au  -  3.424  A 


yielding 


0 Li-Ar  '  3'483  * 


'Ar-Li 


515.66  K 


Since  £2^  changes  very  slightly  for  an  argon  and  lithium  mixture  in 
the  range  of  the  present  calculations,  we  employ  the  average  value 
£2^  =  1.028  for  the  calculations. 

For  the  diffusion  of  lithium  vapor  in  an  argon  layer,  the 


Schmidt  Number  is  given  by 
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